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Abstract. Using a high statistics sample of muon neutrino charged current quasielastic (CCQE) events, we report the 
first measurement of the double differential cross section ( witcos fy, ) ^ or tms P rocess - The result features reduced model 
dependence and supplies the most complete information on neutrino CCQE scattering to date. Measurements of the absolute 
cross section as a function of neutrino energy 

(cT[£ir ]) and the 

single differential cross section ( ~? ) are also provided, 

"Qqe 

largely to facilitate comparison with prior measurements. This data is of particular use for understanding the axial-vector form 
factor of the nucleon as well as improving the simulation of low energy neutrino interactions on nuclear targets, which is of 
particular relevance for experiments searching for neutrino oscillations. 
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CCQE EVENT SELECTION IN MINIBOONE 

The MiniBooNE^ detector, a spherical tank filled with mineral oil, is surrounded by 1280 8" photomultiplier tubes 
(PMTs) to detect Cerenkov light from charged particles 0. In the 19.2/is readout window, a "subevent" is defined as a 
timing cluster of PMT hits. The identification of CCQE interactions relies solely on the detection of the primary 
muon Cerenkov light (first subevent) and the associated decay electron Cerenkov light (second subevent) in these 
events ft 

1: v M +« -► n~ +p 

2: pi- -^e- + v e + V^. K ' 

where each line in this equation identifies the subevent where each process occurs. Therefore, a CCQE candidate is 
characterized with a total of 2 subevents. After cuts, 146070 events are identified from 5.58 x 10 20 protons on target 
collected between August 2002 and December 2005. The cuts are estimated to be 26% efficient at selecting CCQE 
events in a 550 cm radius, with a CCQE purity of 78%. 

The largest background is that from CC single-pion production (CCl7C + ). The CCl7C + interaction, proceeds as, 

1: V n +p{n) -> n~ +p(n) + %+ , 7T+ -^/x+ + v M 

2: [1- -*e~ + v e + V fl (2) 

3: ju+ -^e+ + V e + V M . 

Note this interaction results in total 3 subevents, the primary interaction and 2 muon decays resulting in an electron and 
a positron. Although these events can be removed from the CCQE sample by requiring only one muon decay (a total 
of 2 subevents), there is still a significant number of CC1 7t + events that contribute to the CCQE background because 
one of the muon decays may be missed for various reasons. Among them, n + absorption is a large effect (>40%) with 
large uncertainty (~30%). Additionally, the prediction of CCl7T + backgrounds in the CCQE sample rely on the Rein 
and Sehgal's model [5] and final state interactions (FSIs) in the NUANCE event generator [6] which are not sufficiently 
accurate for a precise background prediction to measure the absolute CCQE cross section. 



1 Ph.D thesis work at Indiana University, Bloomington, IN 

2 The mini-Booster neutrino experiment (MiniBooNE) at Fermi National Accelerator Laboratory (Fermilab) is designed to search for — » V e 
appearance neutrino oscillations (1J|. 

3 Detailed information about the Fermilab Booster neutrino beamline and the MiniBooNE neutrino detector are available elsewhere 1 2, 3]. 




FIGURE 1. (color online). The distribution of events in Qq E for the (a) 2 and (b) 3 subevent samples. Data and MC samples are 
shown along with the individual MC contributions from CCQE, CC1tt + , and other channels. The left side is before the application 
of the CCl?r + background correction. The inset in (b) shows the CC17T+ reweighting function as determined from the background 
fit procedure. The right side is the same distribution after the application of the CC\k + background correction and the new CCQE 
model parameters M e ^ and K as determined from the fit procedure described in the text. 



CC17T+ BACKGROUND MEASUREMENT 

Because of uncertainties in the CCl7T + background predictions, we instead measure the CCl7£ + rate in our CCl7T + 
data and the event generator is adjusted to match. By this, the predicted kinematic distribution of CCl7T + events is 
modified, and the systematic error of CCl7£ + cross section is reduced to the level of the n + absorption uncertainty. 

The left plot in of Figure[T]shows the Qg E distribution^ for data and Monte Carlo (MC) of the two samples before 
the reweighting of CCl7T + MC events. The 2-subevent sample shows good shape agreement between data and MC. 
NUANCE uses the relativistic Fermi gas (RFG) model [7] for CCQE interactions. In the previous work, we adjusted 2 
parameters in RFG model, the effective axial mass M e ^ and Pauli blocking parameter K, to match the shape of the Qg E 
distribution to data |4]. Note that analysis did not consider the overall normalization of events. The 3-subevent sample 
shows a large data-MC disagreement in both shape and normalization. Using these samples, a simultaneous fit was 
performed for the shape and normalization of the 3-subevent sample, and the normalization of the 2-subevent sample. 
These were then used to determine the CC17T+ reweighting function which is shown in the inset plot of FigureQi) (left). 
In order to reduce the sensitivity to the details of the shape of the 2-subevent sample, only the 0.2< Qg E (GeN 2 ) <0.6 

region was considered for the normalization parameter of this function. The Qg E shape of the CCQE sample was 
fit later although it has no impact on the cross section measurements. The effect of the CCQE normalization on the 
3-subevent sample was minimal since the background from CCQE in this Qg E region is small as can be seen in the 
left plot of Figure HI). As a final step, with the measured CCl7T + background incorporated, a shape-only fit to the 
2-subevent (CCQE) sample is performed in order to extract revised CCQE model parameters [4]. The normalization 
of the CCQE sample is then extracted from the fit described above. The Qg E distributions of data from all subevent 

samples is shown together with the MC prediction in the right plot of Figure [T] Data-MC agreement is good in both 

e ff 

subevent samples. A fit to the 2-subevent sample provided adjusted CCQE model parameters, M A and K. This was 
a "shape-only" fit, that is, the MC was normalized with an arbitrary factor to have the same integrated event count as 
the background subtracted data. The fit yielded, 

1.35±0.17GeV/c 2 ; 
1.007 ±0.012; 
47.0/38 . 
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4 The neutrino energy E and 4-momentum transfer Qq e are reconstructed by assuming a CCQE interaction and neutron at rest, with averaged 
nucleon binding energy = 34 MeV. 
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FIGURE 2. (Color online). Left plot is the Q 2 QE distribution of the data, MC before, and MC after the fit with errors. Right plot 

is the I — a contour plot for the — K fit. The filled star shows the best fit point and 1 — <r contour extracted from this work. 
The open star indicates the best fit point and 1 — a contour from the previous work [J]. Two regions are shown from the previous 
work, the larger area indicates the total uncertainty on the results including the background uncertainty @]. 

The left plot of Figure [2] shows the Q 2 q E distribution of data, MC before, and MC after the fit with all sources of 
error. Data and MC after the fit agree within shape errors. The right plot of Fig. [2] is the 1 — a contour regions of 
this fit together with the results from the previous MiniBooNE analysis [4]. Note that the current result is consistent 
(to within 1 — a) with K = 1. This is because the CC17T+ background resulting from the procedure in this work has 
changed by an amount only just consistent with the error assigned on the background in the previous work. The value 
for K is quite sensitive to the CC 1 n + background at low Q 2 q E . However, the previous and current results are consistent 
at the 1 — a level. 

The effect of the new is clearly seen in 2-dimensional plots. Figure [3] shows the data-MC ratio of CCQE 
candidate events as a function of muon kinetic energy 7^ (GeV) and muon scattering angle cos p. Note the muon energy 
and muon scattering angle observables are the basis of all reconstructed kinematics variables in the CCQE channel 
in MiniBooNE. In the left plot, we use the world averaged nuclear parameters (Ml = 1.03 GeV/c 2 , k = 1.000) Jst]. 
As can be seen, data-MC disagreement follows auxiliary lines of equal Q 2 . This is the same tendency observed in the 
previous CCQE analysis in MiniBooNE |4], indicating that data-MC disagreement is more likely due to an incorrect 
cross section prediction (=function of Q 2 ) than an incorrect flux prediction (=function of neutrino energy). After 
introducing the new M4 and K (M^ = 1.35 GeV/c 2 , K = 1.007), Fig. 0right plot, data-MC disagreement is reduced 
and we obtain an improved cross section prediction across the entire kinematic space. 

Note, this modification of the CCQE cross section prediction does not affect the CCQE absolute cross section 
measurement, presented below. 

CCQE ABSOLUTE CROSS SECTION MEASUREMENTS 
Flux-averaged double differential cross section 

Figure. [4] shows the flux-averaged double differential cross section, dT d ~ d ° , for the CCQE process. The flux- 
averaged total cross section, an integral of the double differential cross section (— 1 < cos^ < +1 and < 7^ (GeV) < 
°°) is 9.412 x 10~ 39 cm 2 . The total normalization error on this measurement is 10.8%. 

The kinematic quantities, 7^ and cos^, have been corrected for detector resolution effects only. This result is the 
most model-independent measurement of this process possible with the MiniBooNE detector. No cuts on the recoil 
nucleons are used to define this process. The neutrino flux is an absolute prediction and was not adjusted based on 
measured processes in the MiniBooNE detector. 
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FIGURE 3. Ratio of MiniBooNE CCQE data/simulation as a function of measured muon angle and kinetic energy. Left plot, 



with world averaged l/ff' (=1.03 GeV/c z ) and K (=1.000), and right plot, with newly determined M^ J (=1.35 GeV/c z ) and 
K (=1.007). The ratio forms a 2D surface whose values are represented by the gray scale, shown on the right. If the simulation 
modeled the data perfectly, the ratio would be unity everywhere. Contours of constant E v and Q 2 are overlaid. 
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FIGURE 4. (Color online). The flux-averaged double differential per nucleon (») cross section for the CCQE process. The 
dark bars indicate the measured values and the surrounding lighter bands show the shape error. The overall normalization (scale) 
error is 10.8%. 
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FIGURE 5. (Color online). The flux-averaged single differential per nucleon (n) cross section for the CCQE process. The 
measured values are shown as points with the shape error as shaded bars. Predictions from the NUANCE RFG model with different 
values for the model parameters are shown as histograms. 



Flux-averaged differential cross section 

Figure. [5] shows the flux-averaged single differential cross section, jfS . The reconstructed 4-momentum transfer 

£)g £ depends only upon the (unfolded) quantities and cos fl . 

In addition to the experimental result, Figure [5] also shows the prediction for the CCQE process from the NUANCE 
simulation with three different variations of parameters in the underlying RFG model. The predictions are flux- 
averaged and absolutely normalized. The RFG model is plotted with both the world-averaged CCQE parameters 
(M A = 1 .03 GeV, if = 1 .000) and with the CCQE parameters extracted from this analysis (M A = 1 .35 GeV, K = 1 .007). 
The model using the world-averaged CCQE parameters underpredicts the measured values significantly (by w 30%). 
The model using the CCQE parameters extracted from the shape fit to the MiniBooNE CCQE data are within w 10% 
of the data, consistent within the normalization uncertainty of w 10%. The prediction with the CCQE parameters from 
this analysis scaled by 1 . 10 is also plotted and is in good agreement with the data. 



Flux-unfolded total cross section 

The flux-unfolded total cross section (o[Ey ' ]) as a function of estimated neutrino energy £f ' is shown 

__ OE RFG 

in Figure [6] The quantity Ey ' is a model-dependent estimate of the neutrino energy obtained after correcting for 
both detector and nuclear model resolution effects. These results depend on the details of the nuclear model used for 
the calculation. The dependence is only weak in the peak of the flux distribution but becomes strong at E v < 0.5 GeV 
and E v > 1.0 GeV, in the "tails" of the flux distribution. 

In Figure [6] data are compared with the NUANCE implementation of the RFG model with the world averaged 
parameter values (M^ = 1.03 GeV, K — 1.000), and the parameters extracted from this work (M^ = 1.35 GeV, 
K = 1.007). These are absolute predictions from the model — they are not scaled in any way. The measurement is 
^20% higher than the RFG model prediction with world average parameter values at the flux peak (700 — 800 MeV). 
The prediction with the RFG parameter values extracted from the shape-only fit to MiniBooNE CCQE data reproduces 
the data significantly better, to within 1 a for every point over the entire measured energy range. 

Figure|6|b) shows the CCQE results from the LSND |9] and NOMAD 1 10] experiments. It is interesting to note that 
NOMAD results are better described with the world-average AfJ and K values. 
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FIGURE 6. (Color online). The flux-unfolded total per nucleon (n) cross section with total errors and bin widths plotted indicated 
with the data points. In (a) shape errors are shown as shaded boxes. In (b), a larger energy range is shown along with results from the 
LSND [9] and NOMAD 1 10] experiments. Predictions from the NUANCE simulation with two different RFG parameter variations 
are shown in both plots. 



At this time, a solution to this growing mystery is not evident. Although there are tremendous efforts to model 
this process Hill , no models seem to be able to produce the (1) large observed Mjf and (2) large observed total cross 
section, while keeping the "bare" Ma = 1 .03 GeV (the world averaged value). Model-independent cross section results 
from the MINOS near detector fl^tl . running with E v ^3 GeV and near-future experiments such as MINERvA lfl3ll . 
running with 2 < E v < 20 GeV could help shed further light on this subject. 
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